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The nucleus 54 Mn has been observed in cosmic rays. In astrophysical environments it is fully 
stripped of its atomic electrons and its decay is dominated by the (3~ branch to the 54 Fe ground 
state. Application of 54 Mn based chronometer to study the confinement of the iron group cosmic rays 
requires knowledge of the corresponding halflife, but its measurement is impossible at the present 
time. However, the branching ratio for the related (3 + decay of 54 Mn was determined recently. We 
use the shell model with only a minimal truncation and calculate both /3 + and j3~ decay rates of 
54 Mn. Good agreement for the /3 + branch suggests that the calculated partial halflife of the f3~ 
decay, (4.94 ± 0.06) x 10 5 years, should be reliable. However, this halflife is noticeably shorter than 
the range 1-2 x 10 6 y indicated by the fit based on the 54 Mn abundance in cosmic rays. We also 
evaluate other known unique second forbidden j3 decays from the nuclear p and sd shells ( 10 Be, 22 Na, 
and two decay branches of 26 Al) and show that the shell model can describe them with reasonable 
accuracy as well. 

Pacs Numbers: 23.40.-s, 23.40. He, 21.60.Cs 
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I. INTRODUCTION 

The nucleus 54 Mn decays in the laboratory dominantly 
by electron capture to the 2 + state in 54 Cr with the 
halflife of 312 days. However, as a component of cos- 
mic rays, 54 Mn will be fully stripped of its atomic elec- 
trons, and this mode of decay is therefore impossible. 
The 54 Mn nuclei were in fact detected in cosmic rays us- 
ing the Ulysses spacecraft They offer an attrac- 
tive possibility to use their measured abundance as a 
chronometer for the iron group nuclei (Sc - Ni) in cos- 
mic rays in analogy to the chronometers based on the 
abundances of other long lived isotopes ( 10 Be, 26 Al, and 
36 CI). With them one can, in turn, determine the mean 
density of interstellar matter, a quantity of considerable 
interest. The use of the long lived nuclei as cosmic ray 
chronometers is reviewed in Ref. ||. The importance of 
54 Mn for the understanding of propagation of the iron 
group nuclei that are products of explosive nuclear burn- 
ing have been stressed in Refs. Hj^)- For this program 
to succeed, however, one must know the halflife of the 
stripped 54 Mn (I n = 3+). The decay scheme of 54 Mn 
is shown in Fig. ||; the dashed lines indicate the decay 
paths of the stripped 54 Mn. 

In two recent difficult and elegant experiments the very 
small branching ratio for the f3 + decay to the ground 
state of 54 Cr has been measured: (2.2 ± 0.9) x 10" 9 @] 
and (1.20±0.26) x 10~ 9 @. By taking the weighed mean 
of these values we extract the averaged branching ratio 
of (1.28 ± 0.25) x 10~ 9 . Combining it with the known 
halflife for 54 Mn of 312.3(4) d @, it corresponds to an 
experimental partial f3 + halflife of (6.7± 1.3) x 10 8 years. 
As explained in |l],^]|7| one expects, however, that the 
decay of the fully stripped 54 Mn will be dominated by 



the at present unobservable j3~ decay to the 54 Fe ground 
state. Previously, the partial f3~ halflife was estimated 
assuming that the (3~ and B + form factors are identi- 
cal. Very recently, in Ref. 0, the ratio of the and 
(3 + form factors was calculated using a very truncated 
shell-model and extending it by comparison with similar 
calculations in the sci-shell. The estimated f3~ halflife 
is (6.3 ±1.3) x 10 5 y 0. In this work we will use the 
state of the art shell model and evaluate not only the 
EC decay rate of the normal 54 Mn, but also both de- 
cay branches of the unique second forbidden transitions 
54 Mn(3+) -► 54 Cr(0+) and 54 Mn(3+) -► 54 Fe(0+). By 
comparing the calculated f3 + decay halflife (or branching 
ratio) to the measured one we hope to judge the reliabil- 
ity of the calculation. We then proceed to calculate the 
halflife of the unknown j3~ decay. 

The decays of stripped 54 Mn are unique second forbid- 
den transitions which depend on a single nuclear form 
factor (matrix element). Halflives of several such decays 
in the p shell ( 10 Be) and sd shell ( 22 Na, and two decay 
branches of 26 Al) are known and have been compared to 
the nuclear shell model predictions in Ref. M. For the 
sd shell nuclei, however, only calculations in a severely 
truncated space were performed in ||. Since that time 
computation techniques and programming skills have im- 
proved considerably. Thus, in order to further test our 
ability to describe this kind of weak decays, we repeat 
the analysis ||, using the exact shell model calculations 
without truncation. At the same time the availability of 
new (and different) experimental data for the 10 Be |l(]] 
and 26 Al (ll]] decays make necessary a new comparison 
between experiment and calculations. 

In order to evaluate the decay rate we use the formu- 
lation of fll2|] . The number of particles with momentum 
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p emitted per unit time is: 



N( Pe )d Pe = ^ P 2 eP lF(Z,W e )C(W e )d Pe 



(1) 



where g is the weak coupling constant, p e and W e are 
electron (or positron) momentum and energy, p v is the 
neutrino (or antineutrino) momentum, and Z is the 
atomic number of the daughter nucleus. All momenta 
and energies are in units where the electron mass is unity. 
For the Fermi function F(Z,W e ) we use the tabulated 
values, and the shape factor C(W e ) for the case of the 
unique second forbidden transitions is of the form 
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The nuclear form factor, in turn, is defined as 
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where i denotes the initial state and / the final one; the 
matrix element is reduced with respect to the spin space 
only (Racah convention G| ) ; ± refers to (i^ decay; t± = 
(r x ±iT y )/2, with t+p = n; g A = -1.2599 ± 0.0025 @; 
and R is the nuclear radius (the final expression for 
C(W e ) is obviously independent of R). 

The functions A2 and A 3 are tabulated in Q. Inte- 
grating the rate formula up to the spectrum endpoint we 
obtain the expression for 1/t and, respectively, for the 
halflife (T1/2 = ln(2)r) in terms of the nuclear form fac- 
tor squared. (For the stripped atoms we correct the end- 
point energy accordingly). For the quantity 27r 3 (ln2)/g 2 
we use the value 6146 ± 6 s (m). Note the usual ft value, 
commonly used to characterize a decay, uses the inte- 
grated phase space factor / of Eq. (1), however without 
the constant <? 2 /(27r 3 15 2 )) and the radius factor R 4 . 



II. SHELL MODEL CALCULATIONS 

In our calculation we consider an inert core of 40 Ca 
with the 14 remaining nucleons distributed throughout 
the p/-shell. We use KB3 |^5[ as the residual interaction 
with the single particle energies taken from the 41 Ca ex- 
perimental spectrum. The Hamiltonian is diagonalized 
with the code ANTOINE |ll| using the Lanczos method. 
It is not yet possible to perform a full p/-shell calculation 
but we can come fairly close. Let's denote by / the f-ji^ 
orbit and by r the rest of the pf shell (p 3 / 2 , P1/2 arid 
/ 5 / 2 ). From the calculations of reference [[17| we know 
that one can get a good approximation to the results 
in the full p/-shell calculation when one considers the 
evolution of a given quantity as the number of particles, 
n, allowed to occupy the r-orbits, increases. We have ex- 
tended the calculations of the previous reference allowing 
up to a maximum of n = 7 particles in the r-orbits. The 
TO-scheme dimensions for this calculation are: 17,136,878 



for 54 Cr; 49,302,582 for 54 Mn; 91,848,462 for 54 Fe. Fig- 
ure H shows the variation of the nuclear form-factor for 
the (3 + and (3~ transitions as a function of the truncation 
level n. We have used the harmonic oscillator wave func- 
tions with b = 1.99 fm. From Figure [| it is clear that our 
calculation has already converged for the n — 5 trunca- 
tion. It is also obvious that it would be inappropriate to 
use only the lowest order corrections (n — 2 for /3 _ and 
n = 3 for (3 + ). 

To see how well the calculated wave function repro- 
duces basic characteristics of the ground state of the 
odd-odd nucleus 54 Mn, note that the electric quadrupole 
moment is calculated to be Q = 34 e fm 2 (with ef- 
fective charges = 1.5 and e„ = 0.5), while the ex- 
perimental value is Q = 33 ± 3 e fm 2 . The magnetic 
moment, calculated with the free nucleon gyromagnetic 
factors is /1 = 2.78 fi]\r, while the experimental value is 
/i = 3.2819 ± 0.0013 /iAr. We have also calculated the 
log ft value for the Gamow- Teller electron capture transi- 
tion to the 2+ state in 54 Cr. The calculated log ft = 6.14, 
where we have used the usual quenching factor of 0.76, is 
in good agreement with experimental value of 6.2. Note, 
however, that quenching of this allowed Gamow- Teller 
matrix element is needed to achieve the agreement with 
the experimental rate. (See [ fl8| and references therein 
for the problem of the GT strength quenching.) 

From Eq. (|) we know that the single particle matrix 
elements needed for the evaluation of the form factor 
involve the expectation value of r 2 between the single- 
particle radial wave functions. In the evaluation of this 
quantity we have followed two approximations. First, we 
consider harmonic oscillator wave functions. In this case 
all matrix elements are proportional to the square of the 
length parameter b. We use the prescription of refer- 
ence |l9| ] to determine b from the experimental charge 
radius (r 2 )^ 2 J2(J of the parent nucleus; this leads to 
b = 1.99 fm in 54 Mn. Second, we consider Woods-Saxon 
radial wave functions. They have been obtained using 
the potential well that includes spin-orbit and Coulomb 
terms [^l). The radial parameter of the well has been 
adjusted to reproduce the experimental charge radius. 
The values of the form factor and halflives obtained us- 
ing both methods are listed in Table |. 

We have also evaluated the other known unique second- 
forbidden beta decays: 10 Be(/3") 10 B, 22 Na(/3+) 22 Ne and 
26 Al(/? + ) 26 Mg. For the sd-shell nuclei we consider an 
inert core of 16 and the sd-shell as the valence space. 
These transitions were previously computed in Ref. (^] 
using a truncated shell-model calculation. In our case, 
without truncating the sd shell space, we use the Wilden- 
thal USD effective interaction |g2[ an d determine the ra- 
dial parameters following the procedure outlined in the 
previous paragraph. In the harmonic oscillator approx- 
imation we use b = 1.78 fm for 22 Na and b — 1.81 fm 
for 26 Al. Table | contains the results of our calculations. 
For the decay of 10 Be we reproduce the results of refer- 
ence (9), however our determined b parameter, 1.75 fm, 
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is slightly larger than the one used before (b — 1.68 fin). 
The new experimental value for the halflife [|o| nicely 
agrees with the computed one. 



III. RESULTS AND DISCUSSION 

The top row of Table | shows that our shell model re- 
sult agrees with the measured halflife of the j3 + decay 
54 Mn — > 54 Cr within errors, without quenching of the 
corresponding form factor. The calculated /3~ halflife, 
(4.94 ± 0.06) x 10 5 y (if we arbitrarily take the average 
value between HO and Wood-Saxon calculations), is no- 
ticeably shorter than the range expected in Refs. [jlj||] 
(1-2 x 10 6 y) based on the experimental abundance of 
54 Mn in cosmic rays and the model of the cosmic ray 
confinement. 

Our calculation suggests that the form factor for the 
j3~ decay, 11.7 fm 2 , is larger than the form factor for 
the (3 + decay, 7.8 fm 2 . We can offer some intuitive, al- 
beit very crude, understanding of this difference. Let 
us treat the three nuclei in the extreme single parti- 
cle model. Also, instead of the actual transitions con- 
necting the odd-odd nucleus 54 Mn with the correspond- 
ing even-even ground state, let us consider the transi- 
tions from the seniority zero even-even nuclei to the odd- 
odd one. The Cr — > Mn transition would then involve 
(7^/7/2) 4 (^3/2) 2 -* (tt/7/2) 5 (^3/2), changing a p 3/2 neu- 
tron into a /V/2 proton. In contrast, the Fe — > Mn would 
involve (7r/ 7/2 ) 6 -> (7r/ 7 / 2 ) 5 (1^3/2), changing a / 7/2 pro- 
ton into a p 3 /2 neutron. Using the above naive assign- 
ments, we are led to the conclusion that the j3~ form 
factor should be about \/3 times larger than the (3 + form 
factor. Even though the detailed shell model results are 
not fully determined by the indicated single particle tran- 
sitions (but they are the largest ones), the overall scaling 
factor emerges. 

Among the sd shell transitions in Table |[ the transi- 
tion to the 1.8 MeV state in 26 Mg agrees perfectly with 
the experiment, while the calculated form factors for the 
other two are somewhat larger, by a factor of about 1.5, 
than the experimental value. For 10 Be decay the calcu- 
lated form factor is also a bit larger. We cannot, there- 
fore, draw any conclusion about the necessity of quench- 
ing in the case of the unique second forbidden transitions. 
While the lighter p and sd shell nuclei contain perhaps 
a hint that quenching is needed, it would obviously spoil 
the agreement in the case of the f} + branch of the 54 Mn 
decay. 

In conclusion: Our shell model calculations reproduce 
the experimental halflives of the unique second forbidden 
beta decays within a factor of less than two. No clear evi- 
dence for the quenching of the corresponding form factors 
emerges. For the stripped 54 Mn decays, the shell model 
describes the /3 + branch within errors. It predicts that 
the form factor for the (j~ decay is larger than the one 
for the j3 + decay. The calculated (3~ halflife (and there- 



fore also the total halflife) is noticeably shorter than the 
range based on the observation of 54 Mn in cosmic rays. 
This conflict, albeit relatively mild, makes attempts to 
determine the branching ratio for the (3~ decay experi- 
mentally even more compelling. 



ACKNOWLEDGMENTS 



P.V. is supported in part by the U. S. Department of 
Energy under grant No. DE-FG03-88ER-40397. G.M.P. 
is supported in part by the DGICyES (Spain). Compu- 
tational resources were provided by the Center for Ad- 
vanced Computational Research at Caltech. 



[1] M. A. DuVernois, Phys. Rev. C 54, R2134 (1996). 
[2] M.A. DuVernois, Astrophys. J. 481, 241 (1997). 
[3] J. A. Simpson, Annu. Rev. Nucl. Part. Sci. 33, 323 
(1983). 

[4] J. E. Grove, B. T. Hayes, R. A. Mewaldt, and W. R. 
Weber, Astrophys. J. 377, 680 (1991) 

[5] R. A. Leske, Astrophys. J. 405, 567 (1993). 

[6] K. Zaerpoor et at, Phys. Rev. Lett. 79, 4306 (1997). 

[7] A.H. Wuosmaa et al, Phys. Rev. Lett. 80, 2085 (1998). 

[8] Huo Junde, Nucl. Data Sheets 68, 887 (1993). Data ex- 
tracted using the NNDC On-Line Data Service from the 
ENSDF database, revision March 6, 1998. 

[9] E. K. Warburton, G. T. Garvey, and I. S. Towner, Ann. 
Phys. (N. Y.) 57, 174 (1970). 
[10] F. Ajzenberg-Selove, Nucl. Phys. A 490, 1 (1988). Data 
extracted using the NNDC On-Line Data Service from 
the ENSDF database, revision March 6, 1998. 
[11] P.M. Endt, Nucl. Phys. A 521, 1 (1990). Data extracted 
using the NNDC On-Line Data Service from the ENSDF 
database, revision March 6, 1998. 
[12] H. Behrens and J. Janecke, Numerical Tables for Beta 
Decay and Electron Capture, Landolt-Bornstein, Vol. 1/4, 
Springer, Berlin, 1969. 
[13] A. R. Edmons, Angular momentum in Quantum Mechan- 
ics (Princeton University Press, Princeton, New Jersey, 
1960). 

[14] I. S. Towner and J. C. Hardy, in Symmetries and funda- 
mental interactions in nuclei, edited by Wick C. Hax- 
ton and Ernest M. Henley (World Scientific, Singapore, 
1995). 

[15] A. Poves and A. P. Zuker, Phys. Rep. 70, 235 (1981). 

[16] E. Caurier, code antoine, Strasbourg, 1989. 

[17] E. Caurier, G. Martmez-Pinedo, A. Poves, and A. P. 

Zuker, Phys. Rev. C 52, R1736 (1995). 
[18] E. Caurier, A. Poves, and A. P. Zucker, Phys. Rev. Lett. 

74, 1517 (1990). 
[19] I. S. Towner, J. C. Hardy and M. Harvey, Nucl. Phys. A 

284, 269 (1977). 



3 



[20] H. de Vries, C.W. de Jager, and C. de Vries, At. Data 

Nucl. Data Tables 36, 495 (1987). 
[21] A. Bohr and B. R. Mottelson, Nuclear Structure, Vol. I 

(Benjamin, New York, 1969). 
[22] B.H. Wildenthal, Prog. Part. Nucl. Phys. 11, 5 (1984). 



0+ 



835 




54 Mn * 
' T|,,= 312d 
Q EC = 1377 keV 
Qp-= 697 keV 



54 Fe 



54, 



Cr 




truncation level n. The form factor has been evaluated using 
harmonic oscillator wave functions with b — 1.99 fm 



FIG. 1. Decay scheme of 54 Mn. 



TABLE I. Form factors and half-lives for the (3 and /3 unique second-forbidden transitions. 







B? A F° 21 (fm 2 ) 






halflifc (years) 






HO 


Woods-Saxon 


Exp. 


HO 


Woods-Saxon 


Exp. 


54 Mn(/3 + ) B4 Cr 


7.82 


7.76 


7.1 ± 0.7 


5.55 x 10 s 


5.64 X 10 s 


(6.7 ± 1.3) x 10 s 


54 Mn(/3-) 54 Fo 


11.7 


11.6 




4.89 x 10 5 


4.98 x 10 5 




22 Na(/3+) 22 Nc 


9.24 


9.78 


6.0 ± 0.8 


2.04 x 10 3 


1.87 x 10 3 


(4.8 ± 1.3) x 10 3 


26 A i (/3+) 26 Mg a 


2.44 


2.78 


2.38 ± 0.05 


8.64 x 10 5 


6.65 x 10 5 


(9.1 ± 0.4) x 10 5 


26 Al(EC) 26 Mg a 


2.44 


2.78 


2.39 ± 0.05 


4.58 x 10 6 


3.52 x 10 6 


(4.8 ± 0.2) x 10 6 


26 Al(EC) 26 Mg b 


12.6 


13.8 


8.8 ± 0.5 


1.43 x 10 7 


9.44 x 10 6 


(2.7 ± 0.3) x 10 7 


10 Bc(/3-) 10 B 


23.1 


23.3 


20.4 ± 0.4 


1.18 x 10 6 


1.16 x 10 6 


(1.51 ± 0.06) x 10 6 



a The first-excited state at 1.809 MeV in 26 Mg 
b The second-excited state at 2.938 MeV in 26 Mg 
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